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Recently, the possibility of atomic quantum simulation of dynamical gauge fields was pointed out. 
In this simulation, a system of cold atoms trapped on each link in an optical lattice is regarded as a 
simulator of a U(l) lattice gauge model under certain conditions. However, to implement local gauge 
invariance, a fine tuning of interaction parameters among atoms is necessary. We consider a general 
cold-atom system with this fine tuning being relaxed and investigate the effect of gauge-symmetry 
breaking. We reveal the fact that a wide variety of cold atoms put on the links of an optical lattice 
is still to be a faithful quantum simulator for a lattice U(l) gauge-Higgs model containing a Higgs 
field sitting on site. Physical meaning of the Higgs phase in various quantum simulators is also 
discussed. 
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In the last decade, the idea to use ultracold atoms in 
an optical lattice (OL) as a simulator for various mod- 
els in quantum physics seems to become more and more 
realistic [H, In particular, one of the interesting pos- 
sibilities is to simulate lattice gauge theories by putting 
several kinds of cold atoms on links of an OL in a certain 
rule [3l-[T]|. Several proposals for pure U(l) lattice gauge 
theories (LGT) [l^l were given in Refs.[l-[1] and later 
the proposal was extended to quantum electrodynamics 
with dynamical fermionic matter 0, [l3| and non- Abelian 
gauge models [11.] . 

Shortly after the introduction by Wilson 13], LGT 
has been studied quite extensively mainly in the context 
of high-energy physics both by analytical methods and 
by Monte Carlo simulations, and its various properties 
have been clarified so far. However, the above mentioned 
approach using cold atoms in an OL provides us with 
another interesting method for studying LGT. As an ex- 
ample of expected results, the authors of Ref.0 refer to 
clarification of dynamics of electric strings in the con- 
finement phase. Furthermore, the quantum simulation 
is not suffered from the sign problem in contrast to the 
conventional Monte Carlo simulation of fermions, and it 
also gives important insight to the real-time evolution of 
quark-gluon plasma emerging from a heavy-ion collision. 

One characteristic point of this cold-atom approach is 
that the relation to the gauge system is established only 
under some specific conditions. For example, in Refs.[6i- 
[Tl| , one needs fine tuning of a set of interaction parame- 
ters; in other words, the local gauge symmetry is explicitly 
lost when these parameters deviate from their optimal 
values. 

The above mentioned point naturally poses us a serious 
and important question on the stability of gauge sym- 
metry, and potential subtlety of results of experiments 
of cold atoms as simulators of LGTs, because generally 
the above conditions are not satisfied exactly or easily 
violated in actual cold atomic systems. In this Letter, 
we address this problem semi-quantitatively and exhibit 



the allowed range of violation of above conditions, such 
as the regime of interaction parameters, within which the 
results can be regarded as its own properties of the LGT. 

Explicitly, we shall argue that a U(l) LGT in a four- 
dimensional (4D) lattice coupled with a Higgs field de- 
fined on sites is the target gauge system of atomic quan- 
tum simulation by cold atoms. We take the London limit 
of the Higgs field and ignore its radial degrees of freedom 
for the reason that will become clear shortly. The Higgs 
couplings may be nearest- neighbor (NN) ones and/or 
next-NN ones. These Higgs couplings, as we shall ex- 
plain, effectively describes the effect of "gauge-symmetry 
breaking" terms of the original cold atoms in a 3D OL. 

Let us start with the path-integral representation of 
the partition function Z of the compact U(l) pure LGT, 
the reference system of the present study: 
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where x — {xi,X2, X'i, X4) is the site index of the 3-|-l=4D 
lattice(a;4 is the imaginary time in the path-integral ap- 
proach) and and v (= 1, 2, 3, 4) are the direction indices 
and we use them also as the unit vectors in the /i-th and 
i/-th direction. The angle variable Ox^i G [0, 2ti) and its 
exponential C/xu are the gauge variable defined on the 
link {x, X -\- jS) [l4]. The bar in Ux^ implies complex con- 
jugate, and C2(= 1/e^) is the inverse self gauge coupling 
constant. The product of four Ux^ is invariant under the 
following local (x-dependent) U(l) gauge transformation, 

Ux^ U'x^ = Vx+iMxpYx, Vx = exp(iA^), (2) 

and so are the field strength Ox^u and the action A. 
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It is known [T^l that the system has a weak first-order 
phase transition at C2 = C2c — 1-0. For C2 < C2c the 
system is in the confinement phase in which the fluctua- 
tions of 9xf_i is strong. For C2 > C2c, the system is in the 
Coulomb phase in which the fluctuations of is small 
and Oxfj. describes almost-free massless particles, which 
corresponds to the photon in the electromagnetism [l^ . 

To obtain the quantum Hamiltonian H for Z, let us 
focus on the space-time plaquette term cos 9xi4: in Z with 
the spatial direction index i(= 1,2,3) and rewrite it as 
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where we used the Villain (periodic Gaussian) approxi- 
mation in the first line and the Poisson's summation for- 
mula in the second line. The term lExi^ iOxi — idrExiOxi 
(t is the imaginary time and / = df /dr) shows that the 
integer-valued field Exi on the spatial link (x, x + i) is 
the conjugate momentum of 9xi- Thus the correspond- 
ing operators at spatial site r = (xi, X2, x^) satisfy the 
canonical commutation relation [Ej.i,9r'i'] — —i5„r5ni. 
Eri is the electric field in the electromagnetism but has 
integer eigenvalues due to the compactness (periodicity) 
of A under Bxfj. — > Oxfj. + Stt. 

The integration over 6x4 can be performed as 

G = / J|rf6i:i,4exp(-i^£;„V^6':j;4) = ]J(5q^,o, 

X x,i X 

Qx = 5ZV,S„, (4) 

i 

where we used Y^x.i Exi^iOxA = - Y.x,i ^i^xi-dxi which 
holds for a lattice with periodic boundary condition. 

One may check that the quantum Hamiltonian H cor- 
respondingto Z is just the one given by Kogut and 
Susskind lldl. 
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The second term corresponds to the magnetic energy 
(V X A)"^ in the continuum [iJl- In fact, by inserting 
the complete sets Ie = Ilr,^J2 E^i M^ri}) {{Eri}\ and 
Ig — Ylr il ^^ri ' I {^r i }) ({6'ri} I J between the short- 
time Boltzmann factors exp(— ArTJ) (At = 13/N), one 
may derive the relation, Z = Tr G exp{—/3H), G = 
rir 0^ = '^iEri [I^l • Qr is the generator of time- 
independent gauge transformation and H respects this 
symmetry as [H,Qr] — 0. Qr = is the Gauss's law to 
be imposed for physical states as a constraint. 

Let us turn to the cold atom studies in Refs.0-[ll|, and 
focus on the quantum simulator using Bose-Einstein con- 
densation (BEG) ^] for concreteness. We write the boson 



y' PriexTp{i9ri), whcrc we 
xfj. in Eq.([T]) because the for- 



operator on the link as ip. 
use the same letter 9ri as 9., 
mer shall be identified as the latter. For a homogeneous 
distribution of atoms, we set pri = pa + fjri and take 77^-1 
as fluctuation operator. Then the Hamiltonian may take 
a form as 



^atom - ^E(E^' 



(6) 



and -ffatom residual interactions. We use the coher- 

ent state \{ipn}) and 1= Ilr.i / dPrid9ri ■ \{i'ri}) {{lpri}\ 

to obtain the path-integral expression of .^atom = 
Tr exp(-/3i/atom) as 

Za.o. = / Wlexp [ - Y: ^..V4^„ - A. H..U^x.)_ 
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The first term in the exponent of R.H.S. of Eq.© comes 

from J^xi'^xi'^^'^xi ^ *Sa:4 

and shows -firi 

is the conjugate momentum of 9ri, so Eri — —"Hri- It 

implies that the first term of Eq.® is (27^)"^ 

With Qx = — ^iVxi, we write this factor in Eq.© as 
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This Gaussian factor shows that the constraint Qx = 
for 7 = is shifted to a Gaussian distribution with 
Qlx ^ 7^/At for 7 7^ 0; 7 is a parameter to measure the 
violation of Gauss's law. 

In Eqs. (|7l8p . if one integrates over rjxi S (—00, 00), one 
obtains a term — (4ArVb)~^(V46'a:i — Vi0a:4)^; a part of 
Gaussian Maxwell term. However, this result should be 
improved to respect the periodicity under 9xi 9xi + 2Tr, 
because 9xi is the phase of the condensate. This Gaussian 
term is to be replaced, e.g., by a periodic Gaussian form 
or the corresponding cosine form cos 9xi4 like Eq. ^ (This 
may be achieved by summing over integer rjxi). 

Then we come to consider the following lattice model 
Z~f with compact gauge variables 9xin 



Zj = [dU] exp Y cos 9x1^ + ^ "^2/^1/ cos 9^ 

X,fl x^iiKu 



■^x^y j 1 



(9) 



with ci4 = j'^/At and 02^4 — (2ArVb)~^. Here, the spa- 
tial plaquette C2y-terms in Eq.® may come from i?atom 
and the Ci^-terms are also added for generality mim- 
icking the LGT [There may be other terms in i/atom! ^^e 
Eq.dni)]. 
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The model is equivalent to another LGT with exact 
U(l) gauge invariance, i.e., U(l) Higgs model with asym- 
metric NN Higgs couplings ci^. Its partition function Zi 
is defined as 

C2fj,iy COS Ux^u- 



(10) 



(fix is the phase of the Higgs field cj)x defined on the site 
X with its radial excitation frozen to unity, i.e., (f>x — 
exp{iipx)- So the first term in Ai is the I-shape hopping 
term Cifj.(/)x+fiUxfj.4'x+C-c. Ai is gauge invariant under a 
simultaneous transformation of Eq.Q and 

<px = e'*^" -Kt>'x= Vx(f>x iVx ^ ip'x ^ ipx+ Ax).{ll) 

In fact, Zy is nothing but the gauge fixed version of 
Zi with the so-called unitary gauge ipx — 0. In short, 
the Higgs field (px represents a fictitious charged matter 
field to describe the violation of chargeless Gauss 's law in 
the ultra-cold atoms, where the general Gauss's law with 
charged field is intact. 

The equivalence Z~^ = Z\ is an example of the following 
general rule: Let us start with a general action ^({C/j;^}), 
which is composed of U(l) gauge bit Uxii = exp(i6'2;^) but 
may not be gauge-invariant under Eq.([2]). Then one may 
replace V^^i 4>x+,iUx,^(l>x in A{{Ux,i}) to obtain a Higgs 
action Aii{{Uxf_,}, {(t>x}) = A{{4'x+f^Uxf,(f>x}), which is 
necessarily gauge-invariant under Eqs. ([2]) and (ITlT) . Then 
^{{Uxfi}) may be viewed as a gauge-fixed version of An 
to the gauge (px — I {^x = 0), and their partition func- 
tions Z and Z are equal, i.e.. 



This is proved by the change of variable in R.H.S. as 
Vx^L -> V'x^ = 4>x+^iUxp.(t>x before 4>x integration and using 
= 1, [dW] = [c?[/](Haar measure), and f[d(j)]l = 1. 
In Fig. [1] we show the phase diagram of Zi in the 
C2 — ci plane obtained by Monte Carlo simulation for 
the case ci^ = ci (See Model I there). It shows that the 
confinement and Coulomb phases, which exist in Z of Eq. 
([1]), survives only up to the phase boundary ci — Cic(c2), 
and beyond it the system goes into a new phase, the 
Higgs phase in which both Oxf^ and ipx are stable [l8l |. 
The expectation that the cold atoms may simulate the 
pure gauge theory is assured in qualitative and global 
sense as long as both systems are in the same phases. 
This happens up to the parameters corresponding to ci < 

Clc(c2). 

The three phases can be characterized by the poten- 
tial energy V{r) between two static charges with oppo- 
site signs and separated by the distance r as V{r) oc 
1/r (Coulomb), exp(— TOr)/r (Higgs), r (confinement). 




FIG. 1: (Color online) Phase diagram of four models (|15|l in 
the C2 — Ci,3 plane for the size 16*. The vertical axis IS ci 
for Model I, C3 for Models L and Ls, ci = C3 for Model ItLs. 
All the four models have the three phases. The number (1, 
2) at each critical point indicates its order of transition. The 
Higgs-confinement lines of Model I terminates at C2 ~ 0.7. 



One may distinguish each phase in the experiments of 
cold atoms by measuring atomic density (See Figl2]). 

The above relation p^ . Zj = Zi, leads us to a very 
interesting interpretation that the cold atomic systems 
proposed in Ref-fij and the other related models 0, Q 
have a potentiality to be a simulator of wider range of 
field theories. Namely, the system of cold atoms with a 
general set of values of parameters can be regarded as a 



Higgs 



Coulomb 



Confinement 



FIG. 2: (Color online) Contour plot of deviation of typical 
atomic density Apr = (X)i '7ri/3)^^^ in the xi — X2 plane at 
a;3 = with the external sources of atoms Apcxt = ±pi put 
on the links emanating from r = r± = (±0.4, 0, 0). The white 
regions have Apr greater than certain value and the darker 
regions have the lower Ap,.- The atomic density on the link 
(r,r -\- i) is given by pri = po + rj^i, and the deviation rjxi is 
calculated by using the electric field Eri{= —rjri) with a pair 
of external sources q = ±1 at r = r±. In the Higgs phase, Ap,, 
decreases rapidly away from the sources. In the confinement 
phase, the deviation propagates from one source to the other 
along a one-dimensional string (electric flux). 
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simulator of the U(l) LGT with the Higgs couphngs. 

In all the approaches to quantum simulator of the 
dynamical gauge field presented so far [6l-[lTI|. Hamilto- 
nian of the original system has the form like ffatom = 

(matter). Then, by considering the limit 7^-0 and mak- 
ing the second-order perturbation theory for small g', one 
obtains a gauge-invariant effective Hamiltonian with the 
plaquette C2 term as in Eqs. (I5I9P [gI-IsI. [To| . However, for 

7^0, the original L-shape terms Ulfir+i,j {& -f^Itom) 
remain in the effective Hamiltonian as a first-order con- 
tribution. For large g' , they shall certainly play an im- 
portant role. These L-shape interactions correspond to 
the following terms in the effective action. 



^5 = ^ CSp.L' 



+ C0s{9x+t,,i, - Ox+v,ti) + cos{0xu + Ox+i.,,,) ■ (13) 



Here, cs^i/ = cs^^i, is the coupling between the NN atoms 
on links and = 0. Then, by the same token, this term 
is viewed as a gauge fixed version of the gauge invariant 
term Ai,{{Uxfi.}, {4>x}) obtained by replacing 9xf, (px + 
Ox^i - fx+ti in As. 

The above consideration for cold atoms motivates us 
to introduce the following general Higgs couplings: 



ZiL 

Ail 
Al 



[d4>][dU]expAiLi{Ux^.}AM), 
Ai + Al, 

Cs^iy COs{ipx+i, + Oxp, - 9x1, - Vx+fi) 

x.p<v 

+ COs{ipx + Ox^j. + 6x+i,,i, - (Px+fi+u) 
+ COs{(px+fi + dx+^i^v ~ dx+iy,fi + ^x+v) 

+ COs{ipx+0xv + dx+„.^_,-ipx+i^+p) ■ (14) 



Each term in A^ may be written as 4'x+^i+i'Ux+p,uUxfj.4'x, 
etc. To study the phase structure of Zjl, we consider the 
following cases for definiteness; 

Model £_i ^4 



I 1 1 

L 

Ls 

ItLs 1 



(15) 



Model ItLs (t is for time and s is for space) is the nearest 
to the systems studied in Refs.jl, 0j [fSj - 

Figure[T]is the phase diagram of four Models in Eq.([T5]) 
in the C2 — ci_3 plane. These four phase diagrams have 
similar structure. There are always three phases; Higgs, 
Coulomb, and confinement phases in the order of increas- 
ing size of fluctuations of gauge field Ox^. The order 



of phase transitions is weak-first order or second-order 
in the confinement- Coulomb transition, second-order in 
the Coulomb-Higgs transition, while it changes in the 
confinement-Higgs transition as first-order, second-order, 
and cross over (for Model I) [20] , as C2 decreases from the 
tricritical point at which the three phases merge. 

For sufficiently large C2, Ox^ is almost frozen up to 
gauge transfirmation, and the system reduces to the XY 
model [1^ . Then, for Model I (px has the NN interaction, 
Cifi4)x+p4'x, while for Models L and Ls the interaction 
is next NN one, c^p^(j)x+ij,+v4'x- These XY models ex- 
hibit second-order transition both in 3D (Model Ls) and 
in 4D (Models I, L) between ferromagnetic and param- 
agnetic phases, which corresponds to the Higgs-Coulomb 
transition in FiglTJ The difference of number of interac- 
tion bonds per site, 4, 12, 24 for Models I, Ls, L, makes 
the region of the ordered phase (Higgs phase) larger in 
this order [21*1. 

It is quite instructive to clarify physical meaning of the 
Higgs phase of the gauge system realized in the atomic 
quantum simulators. In the simulator using bosons 
the Higgs phase of the effective gauge system is noth- 
ing but the EEC state as the phase of the bosons (i.e., 
the gauge boson) is stabilized coherently. Therefore, 
the Higgs-confinement transition corresponds to the EEC 
transition. On the other hand, in the quantum simula- 
tor using the Schwinger boson zJT for the gauge field like 
U„ = Zr+iKi^ = 1 or 2) d, [U, the Higgs phase cor- 
responds to the state in which quantum state at each 
link (r, r -I- z) is given by a coherent superposition of the 
particle- number states like |0)r|l)r +i + |l)r|0)r+i- In the 
double-well potential, this state is realized naturally and 
then the Higgs phase of the gauge system appears easily. 

This way of introducing U(l) variables 0, [Il| reminds 
us an approach of starting with the antiferromagnet with 
s — 1/2 quantunis spin at each site and obtaining the 
CP^-l-U(l) LGT [2^, which has a two-component com- 
plex (CP^) variable at each site describing spins and 
an auxiliary hut dynamical U(l) gauge variables on each 
link. Although the CP"'^-fU(l) model and the present 
U(l) Higgs model are different each other, their global 
phase structures are quite similar (See Fig.l of Ref . [22| ) . 

In summary, Fig[T] predicts global phase structures 
of the effective system of cold atoms trapped on 
the links of an OL that are studied in Refs.fl, Q, and 
other similar systems. From the discussion presented in 
Refs.jl, 0] and the relation between the U(l) gauge 
invariant Higgs theory and its gauge fixed theory {^px = 
0) , it may be rather universal that many systems of mul- 
tiplet ( "quantum spins") of cold atoms put on links have 
their U(l) Higgs LGT counterparts. Such equivalence be- 
tween cold atoms and U(l) Higgs model may be called a 
"quantum spin-Higgs gauge correspondence" . In the next 
stage, more explicit and individual relation between an 
cold atom system and its corresponding LGT (jl4p may 
be of interest, which requires expressions of ci^ etc. in 
terms of parameters of the cold atom system. 
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